Effect of dephasing on mesoscopic conductance fluctuations in quantum dots with 

single channel leads 
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We consider the distribution of conductance fluctuations in disordered quantum dots with single 
channel leads. Using a perturbative diagrammatic approach, valid for continuous level spectra, we 
describe dephasing due to processes within the dot by considering two different contributions to the 
level broadening, thus satisfying particle number conservation. Instead of a completely non-Gaussian 
distribution, which occurs for zero dephasing, we find for strong dephasing that the distribution is 
mainly Gaussian with non-universal variance and non-Gaussian tails. 
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I. INTRODUCTION 



Soon after the theoretical prediction of universal con- 
ductance fluctuations (UCF) in disordered electronic 
samples, 1_4 it was shown that their distribution function 
is mainly Gaussian. a These considerations referred to a 
weakly disordered open sample connected to a reservoir 
by broad external contacts. In this case, after diffusing 



through the sample of size L during time r e: 



L 2 /D 



electrons are inelastically scattered in the reservoir (D 
is the diffusion coefficient). The resulting uncertainty in 
the level position, i.e. level broadening 7, is of the order of 
the Thouless energy E c — h/T eTg , and E c 3> A, where A 
is the mean level spacing. This picture remains almost 
unchanged in the presence aLxelatively strong inelastic 
processes within the sample&Q which lead to dephasing 
(the loss of phase-coherence) at the time scale r^ < T crg : 
the conductance distribution remains almost Gaussian 
although its variance is no longer universal but decreases 
oc (T0/T org ) 2-rf / 2 where d is a spatial dimensionality. 

In contrast, it has been shown more recently that the 
conductance distribution is pronouncedly non-Gaussian 
in almost closed systems like chaotic cavities with 



single channel leadsl 
of Coulomb blockadi 



quantum dots in the regime 
where the distribution of the 
heights of the Coulorub. blockade peaks has been mea- 
sured experimentally,E3 or isolated .-conducting rings 
threaded by an Aharonov-Bohm fluxO Weak transmis- 
sion through the contacts means that the electrons typ- 
ically spend more time in the system than r erg so the 
broadening of energy levels 7 <C E c . This inequality 
corresponds to the ergodic regime which allows the use 
of non rp , ert urbative techniques including random matrix 
theoryl 12 a 13 l and the zero dimensional supersymmetric a 
modelo 

There are two distinct types of behavior within the er- 
godic regime, depending on whether the level broadening 
7 is smaller or larger than the mean level spacing A. In 
the absence of level overlapping, 7 <C A, the conduc- 



tance distribution is clearly non-Gaussian. The nature 
of the distribution changes considerably when dephasing 
is substantial, 7 > A, so that the energy levels ovex=. 
lap. This case allows also a perturbative treatmentE£l 
in the framework of the standard diagrammatic expan- 
sion. A crucial question here is how to describe the 
level broadening in quantum dots with point-like exter- 
nal contacts. The first approach, following Ref. |l6|, was 
to attach an additional voltage probe in which electrons 
lose theit,phase coherence before being re-injected into 
the dot .til Another approach was to include an imagi- 
nary potential in the Hamiltonian thus allowdiis-for the 
possibility that electrons ma**-, leave the dot.Lr&Ej How- 
ever, Brouwer and Beenakkerll3 have recently pointed out 
that these approaches do not describe inelastic processes 
within the dot. While the voltage probe accounts for spa- 
tially localized dephasing only, the imaginary potential 
in the Hamiltonian does not conserve the particle num- 
ber. In order to describe dephasing processes within the 
dot, they have introducedEJj in the framework of random 
matrix theory a new variant of the voltage probe model 
which produces spatially uniform dephasing. Thus it was 
found that the conductance distribution became Gaus- 
sian for strong dephasing, r^ <C fr/A. Note that exper- 
iments measuring conductance nuctuations_in quantum 
dots, for example as a function of shapejij also find a 
Gaussian distribution. 

In this paper we present a microscopic diagram- 
matic calculation of the conductance distribution func- 
tion within a dirty quantum dot with two point con- 
tacts in the presence of strong dephasing processes (r^ <C 
h/A) within the dot. We construct a particle-conserving 
model by clearly separating two different mechanisms of 
level broadening: dephasing processes which contribute 
only to the conductance fluctuations and particle leak- 
age which contributes also to the mean value of con- 
ductance. It enables us to resolve a contradiction-Jpp. 
tween the results of the imaginary potential modcfltrEiJ 
that suffered from the violation of particle conservation 



and the voltage-probe model.E-Diia By calculating the mo- 
ments of the conductance distribution, we will show that 
for t^ <C h/A it is almost Gaussian with the variance 



var G =^( G > 2 . 



(1) 



This expression is valid both in the ergodic zero-mode 
regime, r^ 3> T crg , and in the diffusive regime, r# <§C r erg . 
Here (G) is the mean conductance, and (<p(£) = Q<j,(R=C) 
where the small parameter (,,p(R) is effectively a diffu- 
sion propagator at distance R in the presence of strong 
dephasing, and £ is the elastic mean free path. We will 
show that this parameter is given by 



{g 1 ln(L <i> /R), T<j,<.T em , d=2 
{ir/2g )(t/R), r <C r erg , d = 3 
At0/7t, T0 > T crg , any d 



(2) 



Here L,p = (Dr^) 1 / 2 is the dephasing length. Note that 
although the voltage-probe modelta also predicts a Gaus- 
sian conductance distribution, the width of the distribu- 
tion given by Eq. (|f|) is different from that obtained in 
Rcf. nq. In contrast to a cavity with broad leads pd (which 
corresponds to a multi-channel case so that the fluctua- 
tions remain universal) the variance is non-universal even 
in the diffusive regime, t^ -C T crg . 

The factor C<i>(£) which governs the distribution width 
resembles the standard weak-localization parameter. 
This could be seen from the first line of Eq. (0) where 
go ~ E c /A is a classical dimensionless conductance of 
a sample with broad leads. Nevertheless, the origin of 
this parameter is quite different: the variance is not very 
sensitive to the presence of a magnetic field and thus has 
nothing to do with weak localization effects. When time- 
reversal invariance is absent, the variance (|lj) changes 
only by the standard Dyson's factor of 1/(3, typically for 
mesoscopic effects, where (3—1 for the orthogonal en- 
semble (in the presence of potential scattering only), and 
(3 — 2 for the unitary ensemble (in the presence of a finite 
magnetic field or weak scattering by magnetic impurities 
that breaks time-reversal symmetry). 

It is worth noting that the close-to-Gaussian nature of 
the distribution is established in a way somewhat less 
straightforward than for UCF in samples with broad 
leadsJa Had it been possible to restrict considerations 
to the lowest order of perturbation for the variance and 
higher moments, we would have obtained the expig.acB.tial 
distribution as in the imaginary potential modeljZltrEj al- 
beit with a shifted position of the mean. However, the 
leading contribution to the moments is given by the next 
(two-loop) order of perturbation while contributions from 
higher orders are negligible. 



II. TWO KINDS OF DIFFUSONS 

In considerations of conductance fluctuations in quan- 
tum dots with single-channel leads, the lowest order of 



perturbation does not make a leading contribution due 
to the existence of two different perturbative parameters 
related to qualitatively different contributions to the level 
broadening. Although both of them arise also in consid- 
ering any mesoscopic effect for a diffusive sample with 
broad contacts, only one of them is relevant in that case 
as we explain below. _ 

We adopt the impurity diagram techniques to per- 
form ensemble averaging over disorder with allowance for 
dephasing processes. At the first step, we assume a con- 
ventional model of free electrons in a random Gaussian 
potential V, with (V(r)) = and 



(V(r)V(r')) = 



1 



2irvr)T. 



O^cl 



■S(r-r'), 



(3) 



where vq is the one-electron density of states and t c \ is 
the mean elastic scattering time. A standard perturba- 
tive expansions 2 ] is a loop-expansion with the loops being 
built up of diffusons and Cooperons (if time reversal sym- 
metry is preserved). These are impurity ladders (Fig. 1) 
with a pole at small transferred momenta (diffusons) or 
small total momenta (Cooperons). For an isolated d- 
cubic sample of size L with L>l such a diffusion prop- 
agator (diffuson) atT — 0, i.e. in the absence of inelastic 
processes, is given by 



£>(q;w) = 



1 



1 



27TI/07-2 Dqi 



(4) 



where D = VpT e \/d is the diffusion constant, vpT c \ = t, 
and q = nn/ L with n — (m, . . . , nd) being non-negative 
integers, fn the presence of time reversal symmetry, the 
Cooperon propagator has the same form. 

For weak disorder, the diffusion propagator is propor- 
tional to the Fourier transform of the density-density cor- 
relation function, 



K(r;t) = (p(0,0)p(r,t)) , 



(5) 



which describes the probability of finding at the point r 
at time t a diffusive particle injected at the origin at t = 0. 
Particle conservation implies that J K(r;t)d d r — const 
at any time t > 0. On performing the Fourier transform, 
one then obtains that T>(0;u>) = const/a;. 

This exact identity remains valid also at T =/= when 
one includes a weak interaction leading to dephasing pro- 
cesses in the sample. Therefore, the diffuson associated 
with K(r;t) remains the same as in Eq. ([|): the only 
possible effect of dephasing could be in replacing D with 
D c g (for weak disorder, |-D e ff— D\ <C D so that this effect 
may also be neglected). 




x~ 



FIG. bj. A diffuson ladder. A Cooperon ladder is obtained 
by inverting the direction of one of the arrows. 



Mathematically, particle conservation in the diffuson 
associated with the density-density correlation function 
(which we call the intrinsic diffuson) is due to the Ward 
identity which provides an exact cancellation of the one- 
particle self-energy corrections due to inelastic processes 
by vertex corrections (described graphically as interac- 
tion lines between the two sides of the diffusion ladders). 
However, if one defines the diffuson in terms of the sum 
of graphs describing diffusive propagation in the particle- 
hole channel without energy exchange between the par- 
ticle and hole line, it is no longer proportional to the 
density-density correlation function and it is not con- 
strained by particle conservation.^ The absence of the 
usual cancellation of self-energy corrections by vertex cor- 
rections leads ip_the appearance of the cutoff tT 1 in the 
diffusion pole.sH This is precisely the kind of diffuson 
arising in typical lowest-order conductance-fluctuations 
diagrams 1-5 and describing correlations between differ- 
ent members of the ensemble which makes energy ex- 
change between the opposite sides of the diffusion lad- 
der impossible. We will call this kind of diffuson inter- 
sample. Note that all the Cooperon propagators, either 
intrinsic or inter-sample, have exactly the same cutoff, 
tT 1 , as the inter-sample diffusonsJa 

Therefore, one can allow for dephasing by merely in- 
cluding tT as an infrared cutoff in the denominators 
of inter-sample diffusons, i.e. substituting into Eq. (|4|) 
t7 for — iuj at zero frequency. At T ^ 0, averaging 
over energies leadsH'Q to the saturation of the denomina- 
tors of inter-sample diffusons. A consistent description 
of both thermal smearing and_dephasing requires the use 
of the temperature techniques and the development of 
a microscopic theory of dephasing in the partially open 
dot which would provide a dependence of r^ on T and 
microscopic characteristics of the dot. Such a theory is 
well established for relatively large diffusive-samples with 
broad leadsu or for closed quantum dots.c 5 ! 

However, our aim here is just to describe how dephas- 
ing and thermal smearing influence mesoscopic fluctua- 
tions of conductances of the dot. To this end, all we 
need is to distinguish between the two kinds of diffusons: 
inter-sample diffusons which are saturated at q = lu = 0, 
and intrinsic diffusons which diverge at q — lu — 0. Since 
the divergence of this sort could not be treated within 
the diagrammatic approach, we introduce a weak level 
broadening 7 CSC due to escape from the dot which pro- 
vides a cutoff for the intrinsic diffusons. The cutoff in 
the intrinsic diffuson violates the particle conservation 
law but since we restrict our considerations to the case 
7 CSC <C min(T, h/r^), it is dephasing (or thermal smear- 
ing) within the dot rather than particle escape which 
governs the conductance distribution. Even with this 
cutoff, V^q — u = 0) <C T> csc (q = to = 0), i.e. there is 
a sharp distinction between the intrinsic diffusons, T> esc , 
and the inter-sample diffusons, T>^. For simplicity, we 
shall use a zero-temperature diagrammatic technique. If 
T > 7i/t0, one should substitute T for Ti/t^ and the 



thermal-smearing length Lt = (TiD/T) 1 / 2 for the de- 
phasing length L^ in all final results, including that for 
the variance, Eq. ([!]). 

The existence of the two kinds of diffusons was prac- 
tically irrelevant for the description of mesoscopic con- 
ductance fluctuations in samples with broad leads. 1-5 
Firstly, there is strong leakage through the leads, 7 ~ 
h/T eTg , and one could restrict considerations to the case of 
weak dephasing, r erg <C T$, when both intrinsic and inter- 
sample diffusons have the same cutoff. Moreover, even for 
larger samples or higher temperatures when T crg ~S> tj,, 
the intrinsic diffusons do not appear either in (G) or in 
varG but as a third-order perturbation correction which 
is irrelevant in a weakly disordered metal. In contrast 
to this, in the problem under consideration the intrinsic 
diffusons determine the value of (G) and enter varG in 
the leading order of perturbation. 



III. CONDITIONS FOR PERTURBATIVE 
APPROACH 

Traditionally conductance fluctuations in a sys- 
tem with broad, spatially homogeneous contacts are 
considered 1-5 by means of the Kubo formulae! Al- 
ternative-considerations based on the Landauer-Buttikcr 
formulaEl are more convenient in the presence of tunnel 
barriersES or for a lead geometry which involves spatially 
inhomogeneous currents. In this paper we will usCj-this 
formula in the form first derived by Fisher and LeerJ 



G = |rE( r ,£ + r ,2). 



(6) 



-,L(R) 



where the transmission coefficient T ab is the probabil- 
ity of transmission from the channel labeled by a in the 
left (right) lead to the channel labeled by b in the right 
(left) lead. In the case of single-channel leads (of width 
w ~ Tikp ) the point-to-point conductance, G, is given 
by Eq. (|6|) without the summation sign, the transmis- 
sion coefficiexLts in this equation being related to Green's 
functions byEj 



T L(R) 



aia 2 
(hvo¥ 



[« 



+(-) 



ri,r 2 ;e)g-^>(r 2 ,r 1 ;e) , (7) 



where Q + (Q~) is a retarded (advanced) Green's func- 
tion, T\, v-i are the positions of the point contacts, and 
o.\Cti is the transmission probability through the contacts 
themselves. In the entire energy interval of interest the 
mean density of states vq is a constant and the T L ^ are 
energy independent so we will subsequently drop the e 
label. Thus we imply that energy dependence will also 
be irrelevant for the conductance fluctuations: as usual, 
this is valid when the energy difference between two con- 
ductances, lv, is much smaller than the level broadening 
which we assume to be the case. 

We consider the region of parameters defined by 



A < 7 CSC < min {Ti/t^ , T , h/r eIg } , 



(8) 



Here the level broadening due to escape from the dot is 
represented by 7 0SC « A(ai +a 2 ) +7icak where 7i cak is the 
level broadening due to escape from the dot other than 
through the point contacts, for example due to leakage 
through the dot matrix. 

Our perturbative approach is formally applicable both 
to the non-ergodic diffusive regime, r^ < r org , and to 
the ergodic zero-mode regime, r^ > r org . The inequal- 
ity 7osc <§C min(?i/T0, T) is apCraeial distinction from 
the imaginary potential modeC'BEj where the opposite 
inequality holds. This allows us to separate dephasing 
processes inside the dot from those due to weak particle 
leakage. The latter have been introduced only for conver- 
gence of the perturbative approach which is not valid for 
7esc < A. Therefore, a weak violation of particle conser- 
vation is inevitable in the perturbative approach as will 
be explicitly shown later. However, this weak particle 
non-conservation is not related to the level broadening 
due to dephasing inside the dot and, as we shall see, may 
be totally neglected in the final results. On the contrary, 
in the imaginary potential model both the level broad- 
ening and the particle non-conservation have exactly the 
same source and cannot be separated which makes this 
model unsuitable for considerations of the influence of 
inelastic processes inside the dot on the conductance dis- 
tribution. 

The nonperturbative region 7 CSC < A was the 
area of interest for previous zero mode calculations 
However, the exact zero-mode calculations within the 
nonlinear a model cannot be extended (at least, in a 
straightforward way) to include relaxation processes in- 
side the dot by simply introducing r^. The reason is that 
all intra-sample diffusons must and all intrinsic diffusons 
must not contain a cutoff proportional to tT 1 . These 
two kinds of diffusons which are both elementary exci- 
tations in the nonlinear a model must be distinguished 
automatically. This can be easily done in the perturba- 
tive approach to the a modeH by the introduction of an 
extra set of matrix indices to numerate different conduc- 
tance loops. While this additional dependence does not 
lead to any complication in perturbative or renormaliza- 
tion group calculations, it would make direct zero mode 
calculations rather difficult if possible at all. Obviously, 
one circumvents this complication when dephasing is in- 
cluded phenomenologically, as in the voltage probe model, 
where one can proceed with zero-mode calculations O'Ej 
However, there is a price to pay: as two kinds of diffu- 
sons are not distinguished, this model contains, r^ pro- 
portional corrections to the mean conductancetS even in 
the absence of time-reverse invariance when such correc- 
tions cannot appear in all orders of perturbation as we 
will show below. 

The present model overlaps with the voltage probe 
model in the ergodic zero-mode regime (with strong de- 
phasing): r er g <C 70 <C t csc . Even in this case, the over- 
lapping is restricted to the region 7 CSC > A as the re- 



quirement of convergence of perturbation series does not 
allow us to take the limit 7i ca k — >■ 0. Nevertheless, both 
the models address clearly the same physical situation 
where dephasing occurs inside the dot, and the results of 
the present considerations are in a broad agreement with 
the voltage probe model. Moreover, we shall argue later 
that the perturbative approach might be valid even out- 
side the formal limits of applicability. Note finally that 
since our considerations always refer to strong dephasing 
with a noticeable leakage (7i ca k > A), the results are ap- 
plicable both for statistics of the heights of the Coulomb 
blockade peaks (for 0:1,2 <C 1) and for the conductance 
distribution of a single-channel cavity (for a\ t 2 ^ 1). 
Conductance cumulants (((?")) are given by 



«G»» = (e*/h) n (((T L + T*) n )) 



(9) 



where T^ L,R ' should be expressed in terms of Green's 
functions, Eq. (|j), and an extra prefactor of 2 is due to 
spin. Fisher and LeeEj have shown that the unitarity of 
the S matrix implies T R =r-T L even in the absence of 
time reversal invariance (inEH the proof of this equality 
for a two lead system has been explicitly based on par- 
ticle conservation). In the imaginary potential model of 
Refs. Rflfq the particle number is not conserved and, 
if time reversal symmetry is broken, T R 7^ T L as we 
will show explicitly later. That is why the conductance 
distribution in the absence of time reversal symmetry 
obtainedEj from the microscopic foraiula (0) turns out 
to be different from that obtainedErtl from the reduced 
formula G = (e 2 /h)T L . 



IV. MEAN CONDUCTANCE 



In the model under consideration, the particle num- 
ber is almost conserved: for the sake of convergence we 
have introduced 7 0SC > A which describes particle leak- 
age; however, in the region of interest, 7 CSC -C h/r^, the 
leakage effect is small and we will show that in the lead- 
ing order approximation T = T . On average though, 
(T R \ = (T L ) so that the calculation of the mean con- 
ductance is straightforward. 

The mean conductance (G) is clearly proportional to 
the Fourier transform at zero frequency of the density- 
density correlation function, Eq. (0). For the parametric 
region of Eq. (g) and for point contacts separated by a 
large distance, R = |ri — r 2 | ^> t, (G) is dominated by 
the one-diffuson contribution shown in Fig. 2. 





FIG. bl The diagram for the mean conductance. 



Obviously, this is the intrinsic diffuson, obeying par- 
ticle conservation. At each end of the diffuson there are 
two-sided 'petal' shapes given analytically by 

X2 (r,ri) = (g+(v-ri))(g-(r>-r)) (10) 

The filled circles at the petals correspond to the constants 
0:1,2/(2^^0) where a\^ are the transmission probabilities 



through the contacts (here and elsewhere we put fi to 1 in 
all intermediate expressions). The petals correspond to 
motion at ballistic scales since the average Green's func- 
tions (Q ± (r,r')} (drawn as edges of the petal) decay like 
exp(— |r — r'|/2^). At the diffusive scale, |r — r' 



petals may be approximated by 
where the constant \2 is given b 



> £, the 

(r,r') = X2<5(r-r') 



J 



X2 



(Q+(t - r')> (Q-(r' - r)> d d r = / (g+(p)) (<T (p)) d d p = v* 



d^ 



f 2 4- — 



where £ = p 2 /2m — ep. Then one finds 



(G) = 2 A J ^ X2 (r 1 ,ri)P(r' 1 ,r 2 ;7c S c)x2(r 2 ,r 2 )dV 1 d d r 2 = !L ai a 2 ( csc (R) , A = 



2-kvqt . 



a.\a 2 



2tt (2tw ) 2 



(11) 



(12) 



where R = |n — r 2 \ , and ( CS c{R) is defined by 



A 



C( J R)EE2T 2 P(r 1 ,r 2 ;*7) = -^ 



P iq-(n-r 2 ) 



IT l — ' Do 2 + 7 

q 



(13) 



with 7 = 7 osc . As 7 0SC <C E c in the parametric regime 
@; Ccsc(R) is dominated by the zero-mode contribution 
(q = 0) only. Therefore, 



Ccsc(-R) 



A 



7T7e 



Combining Eqs. (121) and (|14|), one finds 



(G) 



e 2 aio^A 

t/i 7esc 



(14) 



(15) 



As expected in the zero-mode approximation, (G) is in- 
dependent of the separation of the point contacts, the 
dimensionality, and the degree of disorder. The perturba- 
tion expansion for (G) is strictly valid only in the regime 
where Ccsc ^C 1. Without introducing 7i ca k, we would 
have 7osc ~ A(ai + a 2 ) < Aso that Ccsc 3> 1. However, 
substituting this value of £esc into Eq. (TL2p reproduces 
the result of the voltage probe model for the classical 
mean conductanceta which corresponds to the standard 
single-channel conductance. 

This means that all perturbative corrections to (G) ei- 
ther vanish or cancel each other. In the absence of time 
reversal symmetry, higher order corrections could result 
from the expansion in the diffuson loops only. In the 
usual case of broad contacts, one should sum over all q 
in Eq. ( |l3| ) which leads to the parameter, g G 7 1 lnL/£ at 
d = 2. It is well known that all the 'main logarithms', 
g~ n hi n cot, are mutually cancelled due.ta the, renormal- 
izability in all orders of perturbation. OE3a The first 
nonvanishing contribution arises in the second order and 
is proportional to g~ 2 \n.ujT i n aor respondence with the 
renormalization group results. r 2 rr 3 However, in the zero 
mode regime considered here there is complete cancella- 
tion at least up to the fourth order similar to that occur- 
ring in the same regime in calculating the energy-level 



correlation functioned. This is easy to verify as the cal- 
culation is much simpler than the standard one in the 
diffusive regime. One could assume that the cancellation 
of the main logarithms in all higher orders in the diffu- 
sive regime corresponds to a similar cancellation of the 
zero mode contributions. Since in the latter case there 
is no renormalization group backing, such a cancellation 
remains only a plausible hypothesis. Nevertheless, it is 
certain that in the absence of time-reversal invariance 
no correction due to dephasing may appear as they are 
forbidden by particle conservation. When time-reversal 
symmetry is not broken, there exist weak localization 
(Cooperon) corrections in the lowest order which are pro- 
portional to the standard g~ x In r^/r in the non-ergodic 
diffusive regime, r^ <C T crg , or to A/r^ in the ergodic 
zero-mode regime, r^ 3> T crg . In any case, in the region 
of Eq. (||) these corrections are negligible, and the main 
result is given by Eq. (Il5|) . 



V. VARIANCE OF CONDUCTANCE 



The lowest order diagrammatic contribution to the 
variance is given in Fig. 3, for clarity both in the stan- 
dard representation 1_4 and in that with diffusion modes 
separated from ballistic onesc 2 ] as in Fig. 1. Although the 
diagram in Fig. 3a looks precisely like one of the leading 
diagrams 1-5 for the UCF, the result of the calculation 
is absolutely different due to the fact that the vertices of 
both conductance loops correspond to the point contacts. 
This is most clear in the representation of Fig. 3b where 
local petals are connected together by two inter-sample 
diffusons. The petals reduce to the same constant as for 
the mean conductance, Eq. (|l l|) . Thus, one finds by anal- 
ogy with Eq. (n2h that this diagram's contribution to the 
variance is 



2A 2 x 2 2 V 2 ( ri ,r 2 ;i 1(t> ). 





FIG. pi The lowest-order contributions to the variance. 
Wavy lines represent diffuson ladders which correspond to 
the inter-sample diffusons. The relation between a wavy line 
and a ladder is the same as for intrinsic diffusons, Fig. oj. 

where the constant A is the overall factor in the expres- 
sion for T L ' R in terms of Green's function, Eq. (0), this 
factor being explicitly defined in Eq. (|l2|). Substituting 
the values of all constants and taking into account the 
contribution of the equivalent Cooperon diagram (ob- 
tained from that in Fig. 3a by inverting the direction 
of arrows in one of the loops) results in the following 
expression 



«G 2 » 



(i) 



1 




l 2 



-aia2(cf>(R) 



(16) 



Here ((f>(R) i s given by Eq. (|13|) with 7 = tI 1 . Its value 
depends on whether one considers the non-ergodic diffu- 
sive regime, r^ <C r crg , or the ergodic zero-mode regime, 
T0 3> T crg . In the former case, the summation in Eq. fll3| ) 
may be approximated by an integration with cutoff at 
q ~ L7 1 = (Dr^) -1 / 2 , while in the latter case only the 
q = term makes a relevant contribution to the sum. 
This leads to the expression for (^ (R) given in the Intro- 
duction, Eq. (En. The contribution ( Jlq ) to the variance 
differs from (G) by the substitution of Cf,(R) for ( CSC (R). 
As C<t>(R) ^ (csc{R) in the region (||), this diffij«j«»ep 
(which is absent in the imaginary potential modclQiJO) 
is crucial for the description of dephasing within the dot. 
It is important that the diffuson and Cooperon dia- 
grams originate from different terms in Eq. ([|). Indeed, 
a diffuson ladder connects Q + and Q~ lines which have 
the opposite directions, while a Cooperon ladder connects 
Q + and Q~ with the same direction. Therefore, in the 
diffuson diagram in Fig. 3a Q + 's have the same direction 
in both loops. Then it follows from the expression for 
j<{L,R) m terms of Green's functions, Eq. (m), that this 
diagram originates from ((T L ) 2 ) + ((T R ) 2 J. As the di- 
rection of arrows must be inverted in one of the loops 
in order to obtain an equivalent Cooperon diagram, in 
this diagram Q +, s have the opposite directions so that it 
originates from 2(T l T r: ). Considering 



rpL _£ rpR £ Qr ^ = 2 (in the absence of time reversal sym- 
metry) when the Cooperon contribution vanishes. For 
f3 = 1 the two contributions in Eq. JjTgLcancel each other. 
In the imaginary potential modcffErEj the lowest order 
diagrams make the leading contributions to the variance. 
Therefore, in this case the reduced formula Qjoc T L is 
not equivalent to the microscopically derived 2 ^ formula 
(g). Since the T L = T R equality is based on particle 
conservation, its breakdown in the imaginary potential 
model is another clear indication of the absence of parti- 
cle conservation. 

In the present model where the particle number is 'al- 
most' conserved, the breakdown of this equality in the 
lowest order indicates that the higher order diagrams 
must be more relevant. This is indeed the case since the 
condition ^ <C Ccsc suggests that the dominant contribu- 
tion to the variance in the region rtq) is given by diagrams 
which contain the lowest possible number of intrinsic dif- 
fusons (Cesc) and have their conductance loops connected 
by only one inter-sample diffuson (£^). 

Two such diagrams are shown in Fig. 4 in the two 
equivalent representations described above. Directions of 
arrows are chosen in such a way that all propagators in 
both the diagrams are diffusons. Inverting the direction 
of arrows in one of the conductance loops in diagrams (a) 
and (c), one converts inter-sample diffusons into Cooper- 
ons, leaving intrinsic diffusons intact. In the diffuson 
diagram (4a), Q +, s have opposite directions in different 
conductance loops so that it originates from 2(T L T R }. 
On the contrary, Q +, s have the same directions in both 
loops in the diagram (4c) so that this diagranr origi- 
nates from <(T L ) 2 ) + ((T R ) 2 ). Thus, in contrast to the 
lowest order diagrams in Fig. Bl even for (3 = 2 when 
time-reversal symmetry is broken, both 2{T L T R j and 
{{T L ) 2 ) + {{T R ) 2 ) contribute to the conductance. 




((T L - T R f) = 2({T L ) 2 ) - 2(T L T R ) , 



(17) 



where, in this lowest diagrammatic order, the first con- 
tribution is given by the diffuson diagram while the sec- 
ond is given by the Cooperon diagram, one can see that 



FIG. kl The leading diagrammatic contribution to the vari- 
ance. In the diagrams b and d, the inter-sample diffusons are 
represented by wavy lines and the intrinsic diffusons by double 
lines. 



The calculation of these diagrams is described in the 
Appendix. For simplicity, we did not draw single- 
impurity lines in diagrams in Fig. 4c and 4d which 
"dress" in a standard wayE 2 ] the inner square in Fig. 4d. 
Such a dressing results in an exact cancellation of one 
of the inter-sample diffusons (wavy lines) by the square 
which makes the contribution of the diagram 4d equal to 
that of the diagram 4b. In addition to these diagrams, 
one can construct only one more contribution to the vari- 
ance which contains two diffuson loops. This contribu- 
tion describes a weak-localization correction to a diffuson 
and contains three intra-sample ladders thus being small 
even compared to the contribution of Eq. (iq). 

Adding together the contributions of all relevant diffu- 
son and Cooperon diagrams in this order, we find 



e 2 a 2 



Cesc(-R) 



«G 2 )) (2) = 
Comparing Eqs. (|l|) and (|l6|), one has 



cm . 



((G 2 )) 



(i) 



c 2 

~ cm 



(18) 



(19) 



in the parametric region (ph corresponding to strong de- 
phasing since (</> (£) is small while (esc is close to 1 . There- 
fore, in this regime the variance is dominated by the non- 
universal contribution of Eq. (jig). Substituting the ex- 



pression for (G) into Eq. fllq), we obtain the result of 
Eq. (jj]) given in the Introduction. 

Let us stress that the exact equality of the diffuson di- 
agrams (|]a) originating from (T L T r ") and (He) originat- 
ing from ((T fl ) 2 ) ensures that the two terms in Eq. Jl7|) 
cancel each other even in the case (3 = 2 when there 
are no Cooperon contributions. The same is valid for 
the leading contribution to the higher order moments 
of T L — T R . Therefore, keeping only the leading con- 
tribution and neglecting sub-dominant contributions to 
((G n )} I {G) n which are due to a small particle leakage, 
one finds that T L — T R as expected in any particle con- 
serving model. 



VI. HIGHER MOMENTS AND DISTRIBUTION 
OF CONDUCTANCE 

In order to determine the distribution, we need to find 
the leading contributions to the nth cumulant. The low- 
est ordep-diagrams are a generalisation of those for the 
variancetJ, Fig. 3. They have two n-petal daisy ver- 
tices, the petals being connected by n inter-sample diffu- 
sion propagators. Each diagram gives a contribution of 
(n - 1)! (G) n [( <p (R)/( csc (R)] n where a factor of (n - 1)! 
arises upon counting all possible ways of ordering the n 
loops. Diagrams which generalize the leading order con- 
tribution to the variance for n — 3 are shown in Fig. H. 






FIG. |5|. Some of the diagrams making a leading contribution to the third moment of conductance. Diagrams in the second 
row are equivalent to their counterparts in the first row. 



In general, such diagrams are obtained by taking n 
conductance loops, each loop having an intra-sample dif- 
fuson (as for the mean, Fig. 2), and connecting these 
loops with the least possible number of inter-sample dif- 
fusons. Actually, the number of inter-sample diffusons 
may vary, as every inner polygon (like the inner square 
in diagram [|<i for the variance) cancels one of them (see 
Appendix) . It is the number of loops consisting of inter- 
sample diffusons which should be minimal. 

For the third cumulant (Fig. 5.), each diagram has two 
such loops. The lowest order diagram, two three-petal 
daisies connected by three inter-sample diffusons, also 
has two diffuson loops. However, there is no cancellation 
there and this diagram is proportional to (<p(R) 3 . Each 
diagram in Fig. 5 is proportional to ClCOCesc- As we 
are interested in the case of strong dephasing ((</, -C 1) 
and weak leakage £ osc ~ 1, this contribution is dominant. 
Note that Fig. 5 does not include all the diagrams of this 
order, and a numerical factor attached to the leading con- 
tribution to (((?")) is not easy to calculate for n = 3 and 
hardly possible for arbitrary n. However, this is not nec- 
essary for what follows. It is sufficient to know that the 
leading contribution to the nth cumulant is proportional 
t° Cj _ W C" sc - Thus, keeping both the lowest order and 
the leading contribution, we find 



((C)) (n- 1)1 

(G) n " 







UR) 



Cosc(R) 



+ C„ 



M*)] r 



(20) 



where n-dependence of c„ is much slower than factorial 
(although it might include a trivial factor like 2" or so). 
For weak particle escape and strong dephasing, Eq. (||), 
the second term in Eq. ( p0| ) dominates. Comparison of 
this dominant contribution to an appropriate power of 
the variance gives 



((G n f 
(var G) n 



mw 2 « 1 



so that the distribution is mainly Gaussian with non- 
universal variance var G given by Eqs. (Q) and ([18]). Note 
that this is similar to the behavior of the higher order 
cumulantsa in the case of UCF in a sample with broad 
contacts, where such a ratio is proportional to g 
Very high cumulants, with n > g^ , grow much—faster 
which leads io lognormal tails of the distributiorc3 sim- 
ilar to thosea for samples with broad contacts. At weak 
disorder the lognormal tails are practically irrelevant as 
they start at 6G ~ (G) y/gS > (G) where SG = G-{G), 
while typical deviations SG ~ yvarG <§C (G). 

In the presence of some particle leakage from the dot, 
C C sc(R) ZL 1; the exponential tails are more important, 
these tails being governed by first term in Eq. ( p0| ) which 
dominates for larger n. It turns out that the exponential 
tails dominate the Gaussian ones for 



(<5G) 2 >varG%^ 



(21) 
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FIG. M. Solid curves show the conductance distribution for 
f3 — 1 with dephasing rates At,/, = 0.1 and 0.025, calculated 
using the second term in Eq. (EOj) with an arbitrarily chosen 
coefficient c„ — 2 n ~ 1 n. For comparison the dashed curves are 
Gaussian distributions with the same mean and variance. 

This estimation results from a comparison of the ex- 
ponential distribution obtained using the first term in 
Eq. (By) with the Gaussian distribution with the vari- 
ance (lq) obtained by neglecting the higher cumulants. 

For weak particle leakage, the exponential tails start 
at deviations which are much stronger than typical ones. 
With increasing leakage, these tails become important 
at the same values of Ccsc at which the variance is 
mainly contributed by escape processes, Eq. (|19[). In- 
equalities ([19h and (Efi) show that the crossover to the 
regime dominated by particle escape happens not at 
t csc ~ max(T0,r crg ) but at considerably weaker parti- 



y/T^jA ;§> T0 in the ergodic 

V Terg 



g in the 



cle escape: either at r et 

regime (r^ > r crg ), or at r osc ~ y/T crg /A > r< 
diffusive regime (t^ ^C r org)- Note that in Eqs. ( |19[) and 
( pl| ) we have disregarded the difference between (^(R) 
and C0W which is exact in the zero- mode regime and 
logarithmically accurate in the diffusive regime at d = 2. 
In the diffusive regime at d = 3, this difference is sig- 
nificant as may be seen from Eq. (0). Taking this into 
account will introduce an extra small factor of ijR into 
the r.h.s. of the estimations ( JlSJ ) and ( |2l| ) thus suppress- 
ing a relative importance of escape processes. 



In the case of main interest, when particle leakage may 
be neglected, the distribution still deviates from Gaus- 
sian since the cumulants, given by the second term in 
Eq. (|20|), do not vanish although they are small. The 
distribution can be restored from these cumulants in a 
standard way (see, e.g., Ref. ||). Exact values of the coef- 
ficient On do not matter as the most relevant contribution 
to the deviation from a Gaussian shape is given by the 
third cumulant, similarly to the case of thfi_conductance 
fluctuations in samples with broad leadso The result- 
ing distributions are shown for two values of t^/t cik in 
Fig. 6. It is seen that deviations from the Gaussian shape 
are not very strong. 



In conclusion, we have found that the conductance 
distribution for disordered quantum dots with single- 
channel leads becomes mainly Gaussian for strong de- 
phasing, 70 -C h/A, in broad agreement with the voltage 
probe model of Ref. |l8[ Let us stress again that since 
we have assumed non-zero 7i ca k in this paper, there is no 
exact region of overlap between our results and those of 
the voltage probe model. Therefore, a direct comparison 
of the results is not possible. However there remain some 
noticeable differences, first of all in the dependence of the 
variance of the Gaussian distribution on r^. If we stretch 
our results to the limit of applicability putting ( csc = 1 
and making the contacts transparent (ax = U2 = 1), 
we will have varG ~ r^ while the result of Ref. [18| for 
transparent contacts is varG ~ t?. Note that although 
in our approach resummation in all powers of Cesc might 
be necessary for C csc ^ 1, it could not change the de- 
pendence of the result on r^. Similarly, in contrast to 
another result of Ref. O, the dephasing rate tT 1 could 
not appear in the mean conductance in the unitary case, 
(3 = 2. Indeed, in this case Cooperons vanish, and parti- 
cle conservation strictly prohibits the appearance of t~7 
in any power of expansion in diffusons. However, this 
disagreement is minor, and is clearly due to differences 
between the models. Let us stress that in contrast to the 
imaginary potential model, our present microscopic con- 
siderations confirm the prediction of the voltage probe 
model on the Gaussian character of the conductance dis- 
tribution with a non-universal variance depending on the 
dephasing rate. 



The most distinctive feature of the perturbative ap- 
proach used here was taking into account two different 
sources of level broadening, one due to the presence of 
contacts entering the classical part of the mean conduc- 
tance, and the other due to dephasing, appearing only in 
weak localization corrections and fluctuations. It remains 
an open question as to whether a similar introduction of 
two different sources of level broadening could be used 
within a non-perturbative approach such as the super- 
symmetric nonlinear a model, thus enabling a calculation 
for 7i cak = 0. 
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APPENDIX A: CALCULATION OF THE 
LEADING CONTRIBUTION TO THE VARIANCE 

The calculation of the diagrams in Fig. [| is straightfor- 
ward in the representation (6) and (d). The contribution 
of diagram ( b) is given by 



SA 2 \ X 3rX2^(ri,ri; »7*)^(ri,r a ; he 



(Al) 



Here the overall factor of 8 takes into account that 
diagram (46) originates from (T L T R \, and there are 
four ways to connect T R and T L with an inter-sample 
diffuson which is clear from the equivalent representa- 
tion of this diagram, Fig. Ep. The constant X3 is con- 
tributed by a triangle consisting of three Green's func- 
tions. This is a constant since, similarly to the calcula- 
tion of the mean conductance, Eq. (Jl2[) , the triangle may 
be approximated at the diffusive scale, |r — r'| > £, by 
X3(r, r', r") = \3 S(r — r') 8(r — r"). The presence of the 
5 fun ctio ns ensures that the contribution of this diagram, 
Eq. ( |Al| ), reduces to a mere product of the local contri- 
butions, X2 and X3, and the diffusons, where the constant 
X3 is given by 



X3= (g + (p)) 2 (g-(p))d d P = -2^ r 2 



(A2) 



At the same scale, the spatial arguments in "D(tx, t\; ij^) 
coincide within an accuracy of order £. Therefore, 
2r 2 D(ri,r 1 ; ij,p) = C, (j) (R = (,) = CM by the definition 



of Eq. Q with (^(R) given by Eq. (g). Substituting all 
the constants into Eq. (Al), we obtain the contribution 
of diagram ( b) as follows 



e 2 aia2 



Cesc(-R) 



W) = CM (Gf 



The contribution of diagram (Ad) is given by 
AA 2 x^V 2 {i lcsc ) J Xi (r) V 2 { ri ,r- i l4> )d d r . 



(A3) 



(A4) 



Here Xi IS contributed by the "dressed" inner square 
which is the sum of one "bare" square and two squares 
containing one single-impurity line e.aGh, this line con- 
necting opposite sides of the square. c3 The factor of 4 
arises because diagram (Ad) originates from ((T L ) 2 ) + 
/ (rpFC^2\j anc j thgj-g are only two ways to connect the con- 
ductance loops with two inter-sample diffusons which is 
clear from the equivalent representation, Fig. He. 



The calculation is more transparent in the ergodic zero- 
mode regime, r org <C t^ <C h/A, when not only the in- 
trinsic diffuson, Eq. (|2j) , but also the inter-sample diffu- 
son, Eq. (|l4|), and the dressed square, X4j reduce to a 
constant. Since the contribution of a single-impurity line 
is given by the correlator (|3j), one finds X4 as follows 



Finally, the overall coefficient for the contribution of 
all the leading diagrams equals 4//3, since in the orthog- 
onal case (P — 1) an equal contribution is made by two 
equivalent Cooperon diagrams. This leads to the result 
ofEq. ©. 



X4= (g + Y(g-yd d p+ 



2lTVQT c 



(Q + Y(Q-)& d p 



The first term is calculated in a similar manner to that in 
Eq. (O) which gives Attuot 3 . The integral in the second 
term coincides with that in Eq. (A2) so that one obtains 



Xi = 47TI/QT 4 



1 \ _ 4tw t 4 

Tel / 1> 



(A5) 



The difference between the one-particle relaxation time, 
r, and the mean elastic scattering time, r e i, is crucial in 
this calculation; the same difference is a formal cause of 
the saturation of the inter-sample diffuson, Eq. (|13|) , for 
q = at 7 = 1/-T0. 

On substituting the zero- mode values of %4, Eq. (|Aq ), 
and 2?(i70), one reduces the integral in Eq. (A4) to 



X4(r)P 2 (r 1 ,r;z7 )d d r=^^ 7 ^ = C0, (A6) 

T <t> \ ZT ) 

so that the dressed square, X4 ; exactly cancels one of the 
inter-sample diffusons. Then, on substituting Eq. ( |A6| ) 
and all the constants into Eq. ( |A4| ) , one finds that in the 
zero-mode case (when (<p(£) = (4, = At^/tt), this con- 
tribution equals that of diagram (46), Eq. (A3). In Eq. 
(A5) we have neglected the 1/t csc contribution to the 
one-particle relaxation time, f/r. Taken into account, 
this would lead to an additional contribution of diagram 
(d) proportional to C^Ccsc rather than C^Ccsc as m -^Q- 
( |A3| ). In this contribution, the inner square cancels one 
of the intrinsic diffusons. Overall, such a contribution of 
diagram (d) would be exactly equal to the contribution 
of the diagram similar to that in Fig. Mb) but made from 
one intrinsic and two inter-sample diffusons. We neglect 
these contributions which are small even compared to 
those of the diagrams in Fig. (g). 

We do not present here the calculation of Xi{ r ) m a 
general case, without the restriction to the zero-mode 
regime. The calculation is standardE3 and results in 



X 4 (r) = 27r^ T 4 [-DV 



l<t> 



so that Xi( r ) is proportional to the diffusion operator. 
Then using the diffusion equation one finds that 



X4(r)T> 2 (r 1 ,r;ij (j) )d d r = 2T 2 V(r 1 ,r 1 ;ij (t ,) 



which generalizes Eq. ( JAq ) and shows that the contribu- 
tion of diagram (MM) in thi s ca se is also exactly equal to 
that of diagram (0») , Eq. ( A3 ) . 
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